The Photon Dispersion as an 
Indicator for New Physics ? 

Wolfgang Bietenholz 

Instituto de Ciencias Nucleares, Universidad Nacional Autdnoma de Mexico, 
A.P. 70-543, CP. 04510 Distrito Federal, Mexico 

Abstract. We first comment on the search for a deviation from the linear photon dispersion relation, 
in particular based on cosmic photons from Gamma Ray Bursts. Then we consider the non- 
commutative space as a theoretical concept that could lead to such a deviation, which would be 
a manifestation of Lorentz Invariance Violation. In particular we review a numerical study of pure 
U(l) gauge theory in a 4d non-commutative space. Starting from a finite lattice, we explore the 
phase diagram and the extrapolation to the continuum and infinite volume. These simultaneous 
limits — taken at fixed non-commutativity — lead to a phase of broken Poincare symmetry, where 
the photon appears to be IR stable, despite a negative IR divergence to one loop. 
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1. LORENTZ INVARIANCE 

Lorentz Invariance is a central concept of relativity: it is a global symmetry in Special 
Relativity, and a local symmetry in General Relativity. Here we stay within the frame- 
work of particle physics as described by quantum field theory, so we refer to Special 
Relativity. Then this symmetry implies that some field <J> (scalar, spinor, vector or ten- 
sor) transforms globally in some representation D of the Lorentz group 50(1,3), 

4>0) ^ C/(A) t 4>(jc)f/(A) = D(A)<J>(A _1 ;c) (AG 50(1,3) , U unitary) . (1) 

A number of high precision tests of Lorentz Invariance involve cosmic rays, see e.g. 
Refs. [1] for recent reviews. They provide perhaps the best chance to probe parameters 
not that far below the Planck scale (in specific cases even exceeding it). Cosmic rays 
attain by far the highest particle energies in the Universe (up to w 10 20 eV). In addition, 
the fact that they travel over tremendous distances may be conclusive for high precision 
properties even at moderate or low energies, because tiny effects could be accumulated 
over a very long trajectory. In view of the latter scenario, we discuss here the photon 
dispersion relation, as a direct kinematic test of Lorentz Invariance. 

2. COSMIC y-RAYS 



Gamma Ray Bursts (GBRs) are emitted in powerful energy eruptions over short periods, 
typically a few seconds or minutes. Temporarily this causes the brightest spots in the 
sky. They were discovered from satellites since the 70s, but nowadays they can also be 



observed from ground. Their sources are small, and there are speculations that they could 
be generated by the merger of neutron stars, or of black holes. Here we are pragmatic 
about their origin: in any case GRBs exist, and we want to see what we can learn from 
the long journey of their photons through the Universe. 

The distance to the source is usually well evaluated from the redshift (if we assume 
the Hubble parameter Hq to be known). In particular, in the year 2005 a spectacular 
GRB was observed [2]: it occurred at a distance of 3.9 Gpc, i.e. it was emitted when the 
Universe was only « 8 • 10 8 years old. Typically the photon energies are in the range of 
10 4 . . . 10 8 eV. Thus it is an obvious idea to use GRBs to probe if the speed of photons 
is really energy independent [3]. 

A systematic study of this question was performed in Ref. [4]. It probed the effective 
dispersion ansatz 

,W(l + £), v(E) = l-£ + ,(£), (2) 

where M is a very heavy mass parameter, which might emerge somehow from some kind 
of "quantum gravity foam". Here we set the speed of light at photon energy E <C \M\ to 
c = 1. If a finite parameter M exists, deviations of the group velocity v from that speed 
should be manifest at very high energy E, or after a very long path. Ref. [4] analyzed data 
of 35 GRBs, detected by three satellites. In fact, photons of higher energy tend to arrive 
slightly later. This could be explained by a relative delay at the source, A SOU rce, which 
would then be amplified by the redshift z to the observed delay A obs = A source ( 1 + z) ■ So 
far this is standard physics. If one includes a finite parameter M, this formula is modified 
to ^ 

— — = A S ource + OuyK(z) , auY = T77T~ , (3) 
l+Z MHq 

for photon energies that differ by AE. Here auv is the Lorentz Invariance Violation 
(LIV) parameter, and the function K(z) is specified in Ref. [4]. The question is now 
whether or not the data for A b s /(1 +z) are compatible with a constant in z (or K(z)). 
The authors of Ref. [4] enhanced the errors by hand until the data became compatible 
with a linear fit. The resulting fit is in fact almost constant, and it suggests 

\M\ > 1.4 • 10 25 eV « 0.001 M Planck (with 95 % C.L.) . (4) 

Some studies of single GRBs or blazar flares even conclude \M\ > 0.01 Mp^n^ [1]. All 
these results are negative regarding the hope to discover new physics, but it is impressive 
that phenomenological information about this energy regime is accessible at all. 

Further hypothetical LIV effects, that are searched for in cosmic photons, are a decay 
7 — > e + + e~, photon splitting, vacuum birefringence and an irregular threshold energy 
for long distance propagation through the cosmic background radiation [1]. 



3. U(l) GAUGE THEORY IN AN NON-COMMUTATIVE SPACE 

Here we address U(l) gauge theory in an non-commutative (NC) space as one specific 
theoretical framework that could lead to a non-linear photon dispersion relation. Due to 
its discontinuous behavior in the IR limit, this approach does no? belong to the large class 



of low energy effective actions — such effective LIV theories, including their predictions 
for photons, have also been reviewed in Refs. [1]. Conceptual issues related to causality 
and energy positivity are discussed in Ref. [5]. 

To describe our framework, let us start by considering a NC Euclidean plane. The 
(Hermitian) position operators obey the commutation relation 

[x^xv] =i® flv = iee flv (ju,v=l,2), (5) 

where we assume the non-commutativity parameter to be constant. This implies a 
spatial uncertainty relation Ax\Ax2 = 0(0). In a Gedankenexperiment, Ref. [6] inter- 
preted this relation (in d = 4) as the event horizon, when one tries to measure extremely 
short distances (of the order of the Planck length), which requires an enormous energy 
concentration in this range. Thus a/II®II could be viewed as a "minimal length", i.e. an 
additional constant of Nature (similar to the parameter 1 /M in Section 2). 

A field theory on such a space is non-local. On the quantum level, the modification of 
the standard geometry at short distances implies not only UV, but also IR effects — in 
particular perturbation theory is in deep trouble due to the emergence of a new type of 
IR divergences ("UV/IR mixing") [7]. 

Here we consider a fully non-perturbative approach, which is even more motivated 
than in commutative quantum field theory. In analogy to that case, we impose a (fuzzy) 
lattice structure of spacing a by means of the operator identity 

/ 2n \ 

exp [i — x^J = 1 . (6) 

The momentum components p^ commute, and they obey the usual periodicity over the 
Brillouin zone, which implies 

e ip^ = e i(Pv+2K/a)^ ^ J = e i(p li +2Tt/a)x ile -\p v x v = ... = ft exp ^JL e ^_ p ^ . 

At fixed 6 and a, this means that the momenta are discrete, and the lattice is automat- 
ically periodic. This is of course in contrast to the commutative space. If the lattice is 
N x N periodic, the momentum components take the form 

In It 
Pn = ^r r nn (n^eZ) = -Na 2 . (7) 

aN % 

Therefore the Double Scaling Limit 

{a ->■ and N ->■ <*>} at Na 2 = const. (8) 

leads to a continuous NC plane of an infinite extent. We start from a finite lattice, 
where both the UV and the IR sector are regular, and we simultaneously remove both 
regularizations in a controlled manner. Other procedures could easily lead to — > or 
6 — > oo, which are different commutative limits. The requirement to link the UV and IR 
extrapolations is related to the notorious UV/IR mixing (for a review, see Ref. [8]). 



A Fourier-type analysis shows that we can return to the use of ordinary coordinates 
Xju, if all field multiplications are turned into star products, such as 

0(x) * yf(x) := 00) exp Q ^ e M vX) V0) • (9) 

The ^-commutator [x^Xy]* :=^*JCv — x v *Xp = i0^ v makes this formulation plausible. 
Hence the action of pure U ( 1 ) gauge theory on a NC plane can be written as 

S[A} = )^jd 2 xF llv -kF llv , F llv = d ll A v -d v A ll +ig{A ll *A v -A v *A ll ) . (10) 

This action is *-gauge invariant (though this does not hold for F^ v ). Note that non- 
commutativity brings in a Yang-Mills type self-interaction term. 

At first sight, it looks hopeless to put this action on a lattice and simulated it: this 
seems to require ★-unitary link variables, Ux u*U Xjfl = 11 (no sum), which can only be 
constructed and updated with stringent constraints given by the complete configuration. 
There is a way out, however, namely the mapping onto a matrix model. 

A long time ago, and in a different context, Gonzalez- Arroyo and Okawa introduced 
the twisted Eguchi-Kawai (TEK) model [9]. It is formulated on a single point, with the 
action 

Stek[U] = Nfi £ Z^Tr[uZulUyU^} (j8 = \/g 2 ) , (11) 

where (jU = 1 , 2 in d = 2) are unitary N x N matrices (which emerged from link 
variables after dimensional reduction of U (N) lattice gauge theory), and Z21 = Z\ 2 = 
exp(27rin/./V) is the twist factor (n e Z, n mod 7^ 0). Much later it turned out that 
this model is equivalent to NC U{\) gauge theory on a. N x N lattice: in particular the 
algebras were shown to be identical (Morita equivalence) [10]. 

Referring to the historic interpretation of the matrices as contracted link variables, it 
is obvious to formulate (rectangular) Wilson loops in the TEK model as 

W^ v (IxJ) = ^zftTtfuyuyuiult] . (12) 

Mapping this term back to the NC gauge theory leads to a sensible definition of a NC 
Wilson loop [11]. (There are no "closed loops" on a NC plane, but the essential criterion 
is *-gauge invariance). A specific Wilson loop is complex in general, but the property 
W^v = Wy^ guarantees that the action is real positive, since it involves a sum over both 
plaquette orientations. Therefore the TEK model formulation is suitable for Monte Carlo 
simulations, i.e. for a non-perturbative study of NC gauge theory. 



4. THE FATE OF THE NON- COMMUTATIVE PHOTON 

Let us now address the question how a possible non-commutativity could deform the 
photon dispersion relation. A 1-loop result of the form 

E2 = P 2 + Jp%2 (C = const.) (13) 



was first derived in Ref. [12]. Hence one could even be tempted to deduce a lower 
bound on ||0|| from the observed high-precision linear dispersion relation (z/0 7^ 0) 
[13]. However, further calculations revealed that the constant C is actually negative, C = 
—2/n 2 [14]. This apparent IR instability is worrisome indeed: it seemed questionable if 
NC QED does have a ground state, i.e. if it actually represents a well-defined quantum 
theory. Hence people quickly proceeded to a supersymmetric version of this model, 
where the negative IR divergence is cancelled [14]. (It is not visible either, of course, 
when one treats the model classically [15].) 

However, we were not happy with this combination of daring hypotheses, so we did 
not assume supersymmetry in addition. We wanted to verify if the NC photon as such is 
really ill-defined [16]. To this end, we considered a "minimally" NC Euclidean space- 
time, which involves a NC plane, [ii ,^2] =i0 = const., plus a commutative plane (x^t) 
which includes the Euclidean time t. (A NC time coordinate would break reflection 
positivity, so in that case the transition to the Minkowski signature would no be on safe 
grounds [17]. Commutative time also alleviates the problems with causality.) 

For the commutative plane we used the standard regularization of an L x L lattice, 
while the NC plane was treated with a TEK model as described in Section 3, i.e. we put 
a N x N matrix model on each of the L 2 lattice sites (L « TV). The first challenge for the 
numerical measurements was the identification of a scale, in order to define the Double 
Scaling Limit (8). A simple ansatz for the lattice spacing, a« 1//3, was successful. It 
corresponds to a Double Scaling Limit which keeps the ratio N/j5 2 constant. Fig. 1 
illustrates, for the example N/[5 2 = 20, that the Wilson loops (W^ v ) of different sizes, 
in various planes, are in fact stable as we vary /V and /3 accordingly. 
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FIGURE 1. Thereat part of square-shaped Wilson loops at a constant ratio N/fi 2 = 20. We show results 
in the commutative plane (left), a mixed plane (center) and the NC plane (right). In the former two planes 
the Wilson loops are real due to the symmetry in the signs of ^3 and t . In the NC plane (W12) has a complex 
phase which rises linearly in the area for large loops. (This Aharonov-Bohm-type effect was also observed 
in 2d NC U(\) gauge theory [18], though with some shape dependence of the Wilson loop expectation 
values [19].) In all cases we observe an excellent Double Scaling behavior. 

Next we measured the open Wilson lines Pnin) in /i -direction of length n in lattice 
unit. In dimensional units this line follows the vector 



&fivPv = nafl (/i : unit vector in /i direction) . (14) 

This observable is *-gauge invariant, in agreement with the non-locality of the *-gauge 
transformation [11]. Such an open Wilson line carries momentum p given in eq. (14). 
It is therefore a suitable order parameter to detect if translation invariance in the NC 
plane is broken. Fig. 2 shows a symmetric behavior at strong coupling (/3 <0.35) and 



again at weak coupling, where the transition value of /3 depends on N. This agrees 
with the expectation based on strong and weak coupling expansions, but in between we 
discovered a phase of broken symmetry, which was not expected from any analytical 
consideration. For the weak coupling transition, Fig. 2 shows a marked hysteresis (the 
two curves for fixed N refer to increasing and decreasing coupling strength), which is 
characteristic for a first order phase transition. 
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FIGURE 2. The measured vacuum expectation values for the open Wilson lines of length 2 and 4 lattice 
units, at N — 15, 25, 35 (from left to right), and j3 fixed by the Double Scaling condition. We observe 
symmetric phases at strong and at weak coupling, but at moderate coupling we discovered a phase of 
spontaneously broken translation symmetry. The weak coupling transition has a clear hysteresis behavior. 

Fig. 3 shows the corresponding phase diagram. For any N the strong coupling transi- 
tion is at /3 ~ 0.35 (cf. Fig. 2), whereas the weak coupling transition occurs at a value 
of /3 N 2 . This means that the Double Scaling Limits — at fixed 6 No 2 N/fi 2 
— all lead to the broken phase, and not to the weak coupling phase as one might have 
expected. Hence it is the broken phase which is really relevant for the question if ob- 
servables in NC U(l) gauge theory are stable or not. The 1-loop result (13) refers to 
the weak coupling phase; it is therefore not relevant for the NC continuum limit (the 
Double Scaling Limit). In the broken phase we found stability for all observables that 
we measured [16]. This holds in particular for our explicit results for the dispersion re- 
lation of the NC photon. It was measured from the exponential decay of the Wilson line 
correlation function, separated in Euclidean time. We set the momentum components in 
the NC plane to zero, hence we measured the energy E(p = p^) \ Pl = P2 =o ■ 

Fig. 4 (on the left) refers to the (symmetric) weak coupling phase. We see an energy 
dip at low momenta. This feature is fully consistent with the prediction (13) based on 
perturbation theory, which some people denote as "tachyonic behavior". (Of course we 
cannot measure negative energies based on the decay of 2-point functions, but the feature 
observed here matches exactly such an IR instability). 

However, as we pointed out above, the phase which really matters for a continuum 
limit to a NC space of infinite extent is the phase of broken Poincare symmetry. The 
result in that phase is illustrated on the right-hand- side of Fig. 4. It follows the standard 
linear dispersion relation. Hence we do not see any non-commutativity effect for the 
dispersion in the commutative plane, but we do see that the disastrous negative IR 
divergence disappears. Our attempts to measure the dispersion relation also at finite 
momentum components in the NC plane were not successful, since we did not obtain an 
exponential decay of the Wilson line correlation functions. 
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FIGURE 3. The phase diagram for pure U(l) gauge theory in 4 dimensions, involving a (spatial) NC 
plane. At fixed size N, Poincare symmetry is intact at strong and at weak coupling, but it is sponta- 
neously broken in a moderate coupling phase. The Double Scaling Limits (DSL) — taken at fixed non- 
commutativity parameter — lead to the broken phase, which is therefore relevant for the NC continuum 
limit. In that phase, the (negative) IR singularity is not observed, and the system could be sensible. 
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FIGURE 4. The dispersion relation of the NC photon in the symmetric weak coupling phase (on the 
left), and in the phase of broken Poincare symmetry (on the right). We show the energy as a function 
of the commutative momentum component p^, E(p = P3)\p i= p 2 =o ■ Our result in the symmetric phase 
is consistent with the IR instability, which was predicted perturbatively [14]. (In contrast, a positive IR 
divergence has been observed in the NC A0 4 model, where translation symmetry breaks as well [20].) 
The physically relevant result, however, is the dispersion in the broken phase, where the Double Scaling 
Limit ends up, and there we observe an IR stable behavior. 



5. CONCLUSIONS 



Lorentz Invariance Violation (LIV) has been searched for intensively in recent years, but 
so far it has not been observed anywhere. We reviewed a failed attempt which checked 
the photon dispersion relation based on Gamma Rays Bursts. 

As one of the theoretical frameworks which could give rise to LIV we discussed 
quantum field theory — in particular pure U(l) gauge theory — in a non-commutative 
space. Our non-perturbative study [16] was based on Monte Carlo simulations in a space 



of Euclidean signature with one (spatial) non-commutative plane. The Double Scaling 
Limit extrapolates simultaneously to the continuum and to infinite volume at a fixed 
non-commutativity parameter. It leads to a phase of spontaneously broken Poincare 
symmetry, which corresponds to a moderate coupling strength in finite volume. 

In this limit we found evidence for convergent observables. This suggests that this 
model might be renormalizable, despite severe problems in its perturbative expansion 
[14]. In particular the photon dispersion in the Double Scaling Limit appears IR stable, 
though we could only measure it in the commutative plane. If one measures a deviation 
from the linear behavior in an NC direction, the confrontation with GRB data could 
establish a robust bound on ||0|| (the norm of the non-commutativity tensor) in Nature. 
So far we provided evidence that the photon may survive in a NC space, although 
perturbation theory found an negative IR divergence (that result only holds in the weak 
coupling phase, which turned out to be irrelevant for the NC continuum limit). 

As an outlook on the phenomenological side, we mention that the Fermi Gamma-ray 
Space Telescope has been in orbit since June 2008. It monitors hundreds of GRBs per 
year, with sensitivity to photon energies E ~ 8 keV...300 GeV [21]. This could fur- 
ther boost the precision of Lorentz Invariance tests on cosmic photons [22], cf. Section 2. 
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